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ABSTRACT 

The  results  of  an  experimental  research  of  the  dependence  of  photoluminescence  (PL) 
intensity  in  region  about  800  nm  for  silicon  nanoinclusions  (quantum  dots)  obtained  by  Si 
ion  implantation  in  Si02  on  the  dose  of  Si  ions  at  two  temperatures  of  an  annealing  T„n  = 
1000  and  1 100°  C are  presented.  It  is  established  that  in  both  cases  the  dependences  have 
the  shape  of  the  curves  with  a maximum.  For  1 100°  C the  maximum  is  shifted  to  the  lower 
dose.  The  influence  of  an  additional  ion  doping  by  the  phosphorus  on  intensity  of  PL  is 
investigated  depending  on  the  dose  (concentration)  of  P and  the  dose  of  the  silicon  at  Tmn 
= 1000°  C.  It  is  shown,  that  in  all  the  investigated  region  of  P doses,  the  presence  of  P 
enhances  the  PL.  The  degree  of  the  enhancement  increases  with  the  P dose,  but  the  rate  of 
the  intensity  enhancement  goes  down.  With  the  growth  of  Si  dose  at  the  constant  dose  of  P, 
the  degree  of  the  enhancement  decreases.  In  an  approximation  of  an  effective  mass,  the 
energy  spectra  of  a quantum  dot  are  calculated  at  the  presence  of  one  or  several  P atoms 
for  various  their  arrangement. 

INTRODUCTION 

The  system  of  silicon  nanoinclusions  (NI)  in  Si02  (Si02:nc-Si)  attracts  the  great  attention 
in  connection  with  the  prospects  of  its  use  in  the  silicon  optoelectronics.  Ion  implantation 
is  one  of  the  most  common  methods  of  the  Si02:nc-Si  production  [1-5].  The  important 
parameter  is  the  intensity  of  the  photoluminescence  (IPL)  for  the  constant  intensity  of  the 
excitation.  This  is  especially  important  taking  into  account  a small  thickness  of  nc-Si  layer 
for  ordinarily  ion  energies.  Generally,  IPL  depends  on  the  several  factors:  NIs  number, 
their  sizes,  the  availability  and  the  concentration  of  the  centers  of  a nonradiative 
recombination.  These  factors  are  sensitive  to  the  conditions  of  the  ion  implantation  (ion 
energy  and  dose),  the  conditions  of  the  annealing  (temperature  and  time),  and  also  the 
presence  of  the  impurities.  In  [6,7]  was  shown  that  for  a fixed  annealing  temperature  TmT1  = 
1000°  C (the  annealing  time,  talln  = 2 hours),  IPL  can  be  adjusted  by  the  choice  of  a dose. 
At  that,  the  sizes  of  NI  practically  do  not  vary,  and  their  density  grows.  The  IPL  growth 
proceeds  up  to  the  dose  for  which  the  NI  joining  begins  due  to  their  high  concentration. 
Also,  it  was  established  that  IPL  can  be  essentially  increased  by  the  additional  implantation 
of  phosphorus.  The  position  of  the  PL  peak  occurs  practically  independent  on  the 
phosphorus  doses.  In  the  report,  the  results  of  a prolongation  of  these  investigations  are 
presented.  The  dependence  of  IPL  on  the  dose  of  silicon  was  taken  at  the  more  high 
temperature  of  an  annealing  - 1 100”  C.  The  influence  of  the  doping  by  P was  investigated 
depending  on  Si  and  P dose.  The  theoretical  research  of  the  influence  of  impurity  atoms 
contained  in  the  quantum  dots  (QD)  on  the  energy  spectra  and  on  the  energies  of  the 
radiative  interband  transitions  is  continued. 
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EXPERIMENTAL 


The  Si02  layers  with  the  thickness  of  300  - 600  nm  were  grown  up  on  the  silicon  in  a “dry 
- wet  - dry”  cycle.  Ion  implantation  by  Si+  and  P*  was  provided  at  the  energy  of  ions  E = 
150  keV  and  the  ion  current  densities  j < 3 |iA/cm2.  The  postannealing  was  carried  out  in 
the  dry  nitrogen  atmosphere.  The  photoluminescence  (PL)  spectra  were  measured  at  the  Ar 
laser  excitation  (X  = 488  nm)  with  a power  density  about  50  W/cm2. 

RESULTS  AND  DISCUSSION 

For  all  the  annealed  after  the  Si+  implantation  layers,  the  PL  spectra  with  a maximum  close 
to  800  nm  are  observed.  Such  spectrum  is  typical  for  Si  Nls  in  SiOi  1 1 -6).  In  Fig.  1 , the  1PL 
dependencies  on  the  Si  dose  at  the  T,,„  = 1000  and  1 100"  C are  presented.  Earlier,  it  was 
shown  [7]  that  for  T.„„  = 1000"  C,  the  dependence  well  coincides  with  the  calculated  one. 
This  calculation  was  made  on  the  base  of  the  model  of  the  Nl  accumulation  with  the  dose 


Fig.l.  The  silicon  dose  (<t>.si)  dependencies  of  PL  peak  intensity  at  two  annealing 
temperatures.  T,„„:  I - 1 000°  C,  2 - 1 1 00"  C. 

and  Nl  joining  at  higher  doses.  It  is  evident  that  for  1 1(H)"  C the  IPL  dependence  on  the 
dose  has  a similar  shape  - the  curve  with  the  maximum,  but  the  position  of  the  maximum 
is  shifted  to  the  lower  doses.  It  is  naturally  to  explain  by  the  fact  that  the  rate  of  Nl 
nucleation  increases  with  the  annealing  temperature,  therefore  for  the  same  concentration 
of  the  excess  silicon  atoms,  the  Nl  density  is  higher.  However,  the  model  used  earlier  for 
Tjnrl  = I (XX)0  C does  not  allow  to  explain  quantitatively  the  IPL  dependence  on  the  dose  for 
Ta™  = 1 100“  C.  Perhaps,  it  is  necessary  to  take  into  account  other  factors,  such  as  Ostwald 
rippening,  the  role  of  which  would  increase  with  Tilr>n. 

In  Fig.2,  the  dependence  of  F = Ip/Lp  (the  ratio  of  IPL  for  Si02:nc-Si  additionally  doped  by 
phosphorus  and  those  versus  without  phosphorus  doping)  on  phosphorus  dose  is  shown. 
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Fig.2.  The  phosphorus  dose  (<DP)  dependence  of  F = Ip/Iwp  (the  ratio  of  IPL  for  Si02:nc-Si 
additionally  doped  by  phosphorus  and  those  versus  without  phosphorus  doping.  3>Si  = 
1T017  cm’2,  Tann  = 1000°  C. 

With  growth  of  the  dose,  the  F magnitude  increases.  In  [6,7],  the  model  is  suggested 
according  to  which  the  IPL  enhancement  at  the  phosphorus  doping  is  connected  with  the 
action  of  two  mechanisms:  the  passivation  of  the  disrupted  bonds  (hence,  the  decrease  of 
the  nonradiative  recombination  rate)  and  the  appearance  of  additional  electrons  in  the 
conductivity  band  of  QD.  In  our  case,  unlike  [5],  the  decrease  of  the  IPL  at  the  high  P 
concentration  is  not  observed. 

In  Fig.3,  the  F dependence  on  the  Si  dose  at  the  constant  P dose  is  presented.  The  decrease 
of  F with  Si  dose  can  be  caused  by  the  reduction  of  P atoms  number  contained  in  each  NI 
with  the  increase  of  the  amount  of  the  latter.  This  is  possible  if  a redistribution  of  the 
phosphorus  between  Si02  and  NIs  during  the  annealing  takes  place. 


<J>si , (1017  cm'2) 

Fig.3.  The  dependence  of  F = Ip/Iwp  (FPL  ratio)  on  the  Si  dose  (®si)-  ®p  = 1-1016  cm'2,  Tan„  = 
1000°  C. 
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THEORETICAL  ACCOUNTS 


To  prove  theoretically  the  absence  of  the  essential  shifts  of  the  PL  peak  at  the  doping  by 
phosphorus,  the  calculations  of  the  electronic  energy  spectra  and  of  the  optical  gap  of  QD 
in  the  envelope  function  approximation  were  carried  out.  For  solving  the  problem  we  were 
made  the  following  assumptions:  I)  the  quantum  dots  have  the  spherical  shape:  2)  the 
impurity  ions  create  a field  of  a hydrogen-like  type;  3)  spin-orbit  interaction  can  be 
neglected;  4)  the  silicon  quantum  dots  are  placed  in  a dielectric  matrix. 

For  the  definition  of  energy  levels  we  must  solve  the  equation  for  the  envelope  functions 

f„(r) 

HtlvFv(r)=£FM(F),  (l) 

where  E is  the  electron  energy  which  is  counted  off  the  top  of  the  valence  band,  and 
operator  H^v  represents  the  single-electron  Hamiltonian  which  was  written  in  a k ■ p - 
approximation.  Near  the  X-point  in  a conduction  band  the  Hamiltonian  is  represented  by  a 
2x2  matrix: 


6 ^HI,  HI, 


Jt2-3fc‘  +V, 


Ji- 


f ■&%- 

xk  -i — -k, 
’ y m, 


\”h  mi 

ft  6 [m,  m,)\  -)  c 


- Arx  + 


is  the  isotropic  operator  with  the  electronic  effective  mass  HJt.  = 3hI,HI, /(2 HI,  + HI,  ), 


where  HI,  and  HI,  are  the  "transverse”  and  "longitudinal"  mass  respectively.  These 
masses  equal  0. 1 f)H0  and  O.O^Hq.  Value  = 0. 1 4 4- 2k/h  is  determing  the 


distance  in  k -space  between  X-point  and  the  nearest  minimum  of  the  energy  in  Brillouin 
zone.  Cl  = 0.54oim  - the  lattice  constant,  Arx  =1.21  .reV  is  the  difference  between 


the  energies  in  X-point  of  the  conduction  band  and  in  T-point  of  the  valence  band.  Vc 
the  operator  of  Coulomb  interaction  which  can  be  written  as 


eA(£s-£,)$Yhi'r'  / + 1 pi  \ 

esR  §£  Rs  +(/  + l)e,(  P'  ; 


is 


g2(£s-£.i)y  ra  1 + 1 

2esR  hR21  /es+(/  + l)e,i 
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It  consists  of  two  part:  the  direct  Coulomb  interaction  between  electron  and  impurity  ions, 
and  the  interaction  between  electron  and  two  image  fields  - of  ions  and  of  electron  itself. 

h j is  the  radius-vector  of  the  impurity  ion  with  the  number  ] , parameters  £s  and  are 
the  dielectric  constants  of  silicon  and  dielectric  matrix,  R is  the  dot  radius,  — C is  the 
electron  charge,  K is  he  number  of  the  impurity  ions  and  P,  (c  O 4) • ) is  a Legendre 
polynomial. 

Near  the  T-point  in  a valence  band  the  Hamiltonian  is  represented  by  a 3x  3matrix 


'H0h+Vx+Vc 
-AMy 
- Nkxkz 


- Nkxky 

H0h+Vy+Vc 

-Nkykz 


-M>y  ' 

- NkyK 

H0h+Vz+VCy 


Here 


h2P 

2m, , ' 


h2  L-M 
2 m0  3 


and  m,t  is  the  hole  effective  mass  which  equals  3w!q/(L  + 2M) . The  numbers  L, 
M,  N are  the  empirical  parameters  which  equal  6.8,  4.43,  8.6 1 respectively  [8], 
The  Hamiltonian  Hyv  takes  into  account  the  real  silicon  band  structure  possessing  rather 

strong  anisotropy.  However,  in  spite  of  this,  the  anisotropy  of  the  dispersion  law  can  be 
described  by  means  of  a perturbation  theory,  as  well  as  a Coulomb  interaction  with 
impurity  ions. 

Solving  equation  (1)  we  find  the  several  lowest  energy  levels  in  the  conduction  band  and 
several  highest  energy  levels  in  the  valence  band  and  also  their  dependence  on  number  K 
for  different  types  of  an  impurity  distribution.  Then  the  change  of  the  optical  gap  A due  to 
the  presence  of  an  impurity  in  the  quantum  dot  can  be  written  as 

8(k,/i)=  a(k,/z)- A(0),  (2) 

where  vector  h means  the  set  {hyh2i...'hK  ). 

The  results  of  our  calculations  are  presented  in  the  Table.  The  gap  ATC  (0)  corresponds  to 
the  case  when  the  potential  barriers  height  is  assumed  infinite.  In  this  case  the  magnitude 
of  the  optical  gap  without  impurity  has  the  greatest  value.  If  we  have  taken  into  account 
the  finite  height  of  the  barriers  (3.1eV  for  electrons  and  4.4eV  for  holes),  the  gap  will  be 
less  ( A^O)  in  the  Table.l).  The  smallest  values  of  the  gap  are  obtained  when  we  take 
into  account  not  only  the  finite  height  of  the  barriers,  but  also  the  discontinuity  of  the 
effective  mass  on  the  boundary  of  the  quantum  dot.  The  gap  A2(0)  was  calculated  for  the 
case  when  the  electron  effective  mass  in  the  dielectric  matrix  equals  WZq  and  does  not 
coinsides  with  the  masses  1Tle  or  ftl,,  in  the  silicon  quantum  dot.  As  we  can  see,  the 
differences  between  AM(0)  and  A2(0)  are  rather  significant.  Note  that  the  values 
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Ao(0)  for  different  sizes  of  the  quantum  dot  are  in  good  agreement  with  the  experimental 
data  [9,10,11]  and  with  the  theoretical  calculations  of  other  authors  who  used  another 
methods  (such  as  local  density  approximation  [12],  tight-binding  model  [13,14]  or 
pseudopotential  method  [15,16]). 

The  calculations  of  the  optical  gap  shift  8(x,/l)  were  carried  out  in  the  framework  of 
more  simple  model  with  the  infinite  barriers.  Another  models  with  the  finite  height  of  the 
barriers  are  significantly  more  complicated  for  the  calculations,  and  as  the  estimations 
have  shown,  the  results  do  not  differ  essentially.  The  position  of  the  energy  levels  and 
hence  the  value  8 will  be  depend  not  only  on  the  quantity  of  donors  but  also  on  the 
distribution  of  ions  inside  the  quantum  dot.  Since  it  is  impossible  to  obtain  the  analitical 
expression  for  the  shift  8 in  a general  case,  we  consider  a three  kinds  of  impurity 
distribution  which  correspond  to  some  extreme  situations. 

Table.  I . The  magnitudes  of  the  optical  gaps  ( A ) and  its  shifts  (8)  due  to  the  presence 
of  the  impurity  atoms  (eV)  for  different  values  of  radius  (nm). 


R 

A,(0) 

A,«» 

A,(0) 

8, 

s2 

K=  1 

K = 3 

vs 

II 

K=  1 

A"  = 3 

K = 5 

1.5 

2.63 

2.28 

1.90 

-0.04 

-0.12 

-0.21 

-0.08 

-0.30 

-0.55 

2 

1.98 

1.83 

1.64 

-0.03 

-0.09 

-0.15 

-0.06 

-0.23 

-0.42 

2.5 

1.69 

1.61 

1.50 

-0.02 

-0.07 

-0.12 

-0.05 

-0.18 

-0.34 

One  of  these  is  the  case,  when  impurity  ions  distributed  by  the  centrally-simmetric  way  on 
the  boundary  of  the  quantum  dot  ( hj  — R ).  As  the  calculations  have  shown,  in  this  case 

8=0  irrespective  of  the  donors  number  K . In  the  other  case,  when  all  the  impurity 
atoms  are  situated  at  the  centre  of  the  quantum  dot  ( hj  =0,  actually  hj  « R),  the 

shift  8,  which  is  denoted  by  8j  in  the  Table,  is  nonzero.  The  last  case  corresponds  to  the 
maximum  asymmetric  distribution  with  the  maximum  electrical  dipole  momentum,  when 
all  the  impurity  atoms  are  situated  in  one  place  on  the  boundary  of  the  quantum  dot.  This 
distribution  type  gives  more  significant  shift  82  (sec  the  Tabic.  1). 

As  it  has  been  seen  from  the  table  the  magnitude  of  the  shift  is  rather  small  almost  in  all 
cases,  shown  in  the  Table,  excluding  the  case  of  the  maximum  asymmetric  distribution 
with  the  large  number  of  donors  ( K — 5).  However,  it  is  clear,  that  this  kind  of  the 
distribution  is  hardly  realized  in  experiments  and  the  possibility  of  such  configurations 
strongly  decreases  with  increasing  K . Hence,  the  contribution  of  the  asimmetric 
configurations  with  the  maximum  dipole  momentum  will  be  negligible. 

In  reality,  the  impurity  atoms,  perhaps,  are  randomly  situated  inside  the  NI.  Therefore,  the 
true  magnitude  of  8 for  K < 5 woukd  be  somehow  averedged  one  between  ~ — 0.2eV 
and  zero,  e.g.  — 0,  ]eV  < 8 < 0.  For  the  case  K » 1 , the  consideration  becomes  not 
justified  enough  because:  (i)  the  interaction  between  impurity  atoms  must  be  taken  into 
account  (the  approximation  of  isolated  atoms  can  be  not  right):  (ii)  the  validity  of  the 
perturbation  theory  can  be  violated. 

Thus,  the  magnitudeS,  obtained  experimentally,  can  not  exceed  a several  per  cents  from 
the  magnitude  of  the  gap  and  the  shift  of  the  luminescence  peak  in  the  presence  of 
phosphorus  must  be  rather  small. 
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CONCLUSIONS 


1.  The  dependences  of  IPL  from  Si  doses  similar  character  for  Ta„„=1000°C  and  1 100°C. 
However,  for  the  latter  case,  unlike  the  first  one,  it  can  not  be  quantitively  described 
by  the  same  model. 

2.  The  additional  phosphorus  implantation  causes  the  IPL  enhancement  for  all  the 
impurity  doses  investigated. 

3.  The  effectivity  of  phosphorus  related  IPL  enhancement  drops  with  the  silicon  dose 
(excess  silicon  concentration) 

4.  Theoretical  consideration  in  effective  mass  approximation  shows  that  (in  agreement 
with  experiment)  the  phosphorus  doping  would  not  shift  the  luminescent  peak  position 
significantly  for  the  interband  transitions. 
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